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Abstract. We utilize single-crystal neutron diffraction to study the C2/m structure 

^ ' of potassium hydrogen carbonate (KIICO3) and macroscopic quantum entanglement 

2~I ' above the phase transition at Tc — 318 K. Whereas split atom sites could be due to 

jy-N , disorder, the diffraction pattern at 340 K evidences macroscopic proton states identical 

to those previously observed below Tc by F. Fillaux et al., (2006 J. Phys.: Condens. 
[^^ . Matter 18 3229). We propose a theoretical framework for decoherence-free proton 

^^ '_ states and the calculated differential cross-section accords with observations. The 

structural transition occurs from one ordered P2i/a structure (T < Tc) to another 
ordered C2/m structure. There is no breakdown of the quantum regime. It is suggested 
that the crystal is a macroscopic quantum object which can be represented by a state 
S^ ■ vector. Raman spectroscopy and quasi-elastic neutron scattering suggest that the 

H I \C2/m) state vector is a superposition of the state vectors for two P2i/a-like structures 

symmetric with respect to (a, c) planes. 
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1. Introduction 

In principle, a crystal with perfect translational invariance is a macroscopic quantum 
systems with discrete phonon states at any temperature below melting or decomposition. 
However, for crystals containing 0-H- ■ -O hydrogen bonds, there is a widely-spread 
dichotomy of interpretation for the so-called "proton disorder", usually observed as 
a dynamical equilibrium between two sites, say 01 — H • ■ ■ 02 and 01 ■ • • H — 02. 
For such systems, semiclassical approaches are commonly used to rationalize proton 
correlation times measured with solid-state NMR and quasi-elastic neutron scattering 
(QENS) [1]. Semiclassical protons are dimensionless particles, with definite positions 
and momentums, moving in a double-well coupled to an incoherent thermal bath. 
Dynamics are rationalized with uncorrelated jumps over the barrier and "incoherent 
tunneling" through the barrier. In fact, these models describe a liquid-like surroundings 
and strong interaction with the thermal bath is supposed to lead to fast decoherence. By 
contrast, diffraction and vibrational spectra provide unquestionable evidences that the 
translational invariance and the quantum nature of lattice dynamics are not destroyed 
by proton transfer. 

The crystal of potassium hydrogen carbonate, KIICO3, is at variance with this 
dichotomy of interpretation [2l [3l HI [5] . This crystal is composed of centrosymmetric 
hydrogen bonded dimers (HC0j)2, separated by potassium ions. Below ^ 150 K, 
protons occupy sites corresponding to a unique configuration for dimers, say L. Between 
150 and 300 K, protons are progressively transferred to the other centrosymmetric 
configuration i?, up to a population ratio i? : L ^ 18 : 82. In contrast to semiclassical 
models, it has been observed that the sublattice of protons is a macroscopically entangled 
quantum object over the whole temperature range. Quantum correlations outside 
any description offered by classical physics are evidenced with neutrons because the 
intensity diffracted by the sublattice is proportional to the total nuclear cross-section 
aH ~ 82.02 bn, instead of auc ~ 1-76 bn for uncorrelated protons. Enhanced diffraction 
observed from cryogenic to room temperatures is consistent with a superposition of 
decoherence-free macroscopic states at thermal equilibrium and the correlation time 
measured with QENS [6] can be rationalized with quantum interferences [5]. There 
are therefore consistent evidences that protons in the crystal field are not semiclassical 
particles possessing properties on their own right [7]. 

In the present paper we examine whether the macroscopic entanglement is destroyed 
by the reversible structural phase transition at T^ = 318 K [8]. The population ratio 
R : L increases discontinuously to 0.5 : 0.5 and then remains a constant above Tc. This 
transition, supposed to be of the order-disorder type P, HOl HU US, UHl HH |15], [16], fTTj . 
could be of great significance to experimental studies of the quantum-to-classical 
interface in the crystalline state. However, structural and dynamical models are still 
in debate and require further studies. We have therefore conducted extensive single- 
crystal neutron diffraction measurements in order to (i) determine the best structural 
model; (ii) observe whether enhanced features due to quantum correlations survive or 
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not above T^. The experiments presented below show that the structural transition 
does not destroy the macroscopic quantum behaviour of protons. It transpires that the 
transition occurs from one ordered structure to another ordered structure through a 
coherent quantum process preserving the translational invariance of the crystal. 

The organization of this paper is as follows. In Sec. Ej we determine the structure at 
340 K, for the sake of comparison with previous models. In Sec. [3], we examine the total 
differential cross-section measured over a large Q-range, at the same temperature, which 
evidence the enhanced features of quantum correlations. The theoretical framework is 
presented in Sec. H] and we calculate the differential cross-section of correlated protons 
in Sec. |5] for comparison with observations. In Sec. [6], macroscopic states are opposed 
to statistical disorder and we utilize Raman spectroscopy to probe proton dynamics. 
We suggest that the crystal is in a pure state that can be represented by state vectors. 

2. Crystal structure 

The reversible structural transition has a very small enthalpy and the crystal retains 
its optical transparency [8]. Below T^, the crystal symmetry is Pixja with four 
indistinguishable KHCO3 entities in the unit cell [3, HI [HI [I9] . Planar centrosymmetric 
dimers (HC0j)2, with O- ■ -O distances of ~ 2.58 A, are practically parallel to (103) 
planes. They are alternately slightly tilted with respect to (a, c) glide planes. 

Kashida and Yamamoto (KY) [9] have measured X-ray diffraction across the 
transition. The unit cell parameters are very little affected, the symmetry above T^ 
is Cljm and (a, c) glide planes become mirror planes. These authors envisaged two 
structural models: (i) the "displacive" model based on reorientation of dimers whose 
symmetry changes from Ci (T < T^ to C2/1 above Tc\ (ii) the "disorder" model composed 
of split Ci dimers, symmetric with respect to the mirror plane. This latter model was 
preferred by KY, but Machida et al. ^1] preferred C^h dimers to fit their neutron 
diffraction data. 

We have collected single- crystal neutron diffraction data with the four-circle 
diffractometer 5C2 at the Orphee reactor (Laboratoire Leon-Brillouin) [20]. The 
twining-free crystal (^3x3x3 mm^), previously studied in Ref. [1] was loaded into 
an aluminum container heated at (340 ± 1) K. This temperature is sufficiently below 
decomposition (?« 400 K) to avoid damage, and it is sufficiently high to marginalize 
critical phenomena [6l [H [121 [131 [13 • '^^^ data reduction was carried out with 
CRYSTALS [211 [22]. 

We conffim the monoclinic C2/m {C^h) symmetry and previous unit cell parameters 
(Table [1]) |ni[ll]- Both structural models, referred to as "C2h" and "C/' in Table [H were 
refined. C2h dimers are composed of four O1/2C1/2OII1/2 entities per unit cell, while 
Ci dimers are obtained with four (HC03)i/2 units. This latter model was refined in 
two steps. First, the split geometry was constrained to fit the mirror-plane symmetry. 
Then, constraints were removed in order to check the stability. After further refinement, 
deviations were insignificant. The final i?-factor is significantly better for the — therefore 
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preferred — Ci model (Fig. [T]and Table [T]), in accordance with KY, but at variance with 
Machida et al. The i?-factor is slightly improved by further splitting protons in two 
sub-sites on either sides of dimer planes, separated by ~ 0.05 A. However, the large 
overlap of the thermal ellipsoids casts doubt on this splitting, so it is ignored in the 
remainder of this paper. 

Positional and thermal parameters for heavy atoms in Tables |2] and [3] are quite 
consistent, within estimated errors, with those of KY, but more accurate. By contrast, 
the proton positions of KY, x/a = 0.0323(32) and z/c = —0.4102(123), are at variance 
with those in Table [2] and the X-ray based O — H bond length of 0.762(52) A is too 
short. These discrepancies could arise from the difference between the maximum of the 
electronic and nucleus density. On the other hand, O — H bond lengths in Table H] are 
virtually identical to those of Machida et al. for the C2h model, and we have checked 
that the C2/1 and Ci models refined with our data lead to identical proton positions. 
Furthermore, proton positions are identical to those previously reported at 14 K and 
300 K [31 H] , although with different occupancies. 

The phase transition can be featured as follows: (i) the discontinuous change of 
the occupancy of proton sites, with no significant change in positions; (ii) small changes 
of K and O atom positions by ~ 0.1 A parallel to (6); (iii) the splitting of carbonate 
groups. The split carbon atoms separated by ~ 0.2 A along (6) can be distinguished, 
while the split oxygen atoms are practically indiscernible. The CO single and double 
bond lengths (1.395 and 1.210 A, respectively in Table H]) are markedly different from 
one another and from those at 300 K (1.333 and 1.273 A, respectively [3]). 

3. Total neutron scattering: measurements 

In order to observe whether the enhanced features previously attributed to macroscopic 
entanglement survive above T^, we have conducted diffraction experiments with the 
SXD instrument at the ISIS pulsed neutron source [23ll2lj. The twinning-free crystal 
previously measured at 300 K (?« 15 x 10 x 3 mm^) was loaded in a vacuum tank at 
(340 ± 1) K. It was shown that multiple scattering effects are negligible [1]. 

Compared to the four-circle diffractometer, the advantage of SXD is twofold. 
Firstly, with the high flux of epithermal neutrons at the spallation source, one can 
measure a more extended domain of reciprocal space than at a reactor source. Secondly, 
with the time-of-fiight technique, the total scattering function is measured at once and 
signals in addition to Bragg peaks are effectively collected. 

By definition, Q^, Qy, Q^, are projections of Q onto the normal mode axes for 
protons, X, y, z (Fig. [1]). The maps of intensity (Figs [2] and [3]) for {a*^c*^ planes at 
different fc-values ik, = Qy/b*, with b* ~ 1.112 A~^) are concatenations of measurements 
for 8 different crystal orientations, with (6) || (6*) either parallel or perpendicular to the 
equatorial plane of the detector bank. 

The anisotropic incoherent scattering intensity centered at Q = 0, that was clearly 
observed at 15 K [3], is barely visible at 340 K. The shortfall of intensity with |Q| is a 
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consequence, via the Debye- Waller factor, of the thermal population of the density-of- 
states [11[25]. 

Cigar-like shaped ridges of intensity along Qz are visible in Fig. [2] for particular 
A;- values. For fc = 0, they appear at Qx = and ±(10.00 ±0.25) A~^ For A; = 1, they are 
barely visible. For fc = 2 or 3, we observe a different pattern of ridges at Qx = ±(5 ±0.2) 
and ±(15 ±0.2) A^^, still along Q^. These features are best visible in Fig. [3]for k = 2.5. 
There is no visible ridge at fc = 4. Then, from fc = 5 to 9, we observe the same sequence 
as for A; = to 4. These ridges are basically identical to those previously reported at 
300 K or below for the ordered P2i/a crystal [3l H]. Consequently, they are not due to 
transition induced phenomena, such as disorder, or domains, or planar faults, or critical 
behaviour, for which coherent diffuse scattering should vanish below T^ |26]. Below it is 
shown that these features can be attributed to diffraction by the sublattice of entangled 
protons. 

4. Theory 

4-.1. The adiabatic separation 

Within the framework of the Born-Oppenheimer approximation, the vibrational 
Hamiltonian can be partitioned as 

?iv = 'Hu + 'Ha.t + ^Hat, (l) 

where Tin and Tiat represent the sublattices of protons (H+) and heavy atoms, 
respectively, while Cnat couples the subsystems. For OH- ■ -O hydrogen bonds, coupling 
terms between OH and O- ■ -O degrees of freedom are rather large [271 EH], and beyond 
the framework of the perturbation theory. Two approaches, either semiclassical or 
quantum, are commonly envisaged. 

Semiclassical protons are dimensionless particles, with definite positions and 
momentums, moving across a potential hypersurface [291 EOl EH |32l |33]. Complex 
trajectories involving heavy atoms lead to mass renormalization, and to incoherent 
phonon-assisted tunnelling [H [6l [31] . This approach is quite natural when the Born- 
Oppenheimer surface is calculated from first principles, but quantum effects are severely 
underestimated. 

Alternatively, if the classical concept of "trajectory" is abandoned, adiabatic 
separation of Tin and 7ia.t may lead to tractable models [H EH EH ES, ESI EZj. Then, 
protons in a definite eigen state should remain in the same state in the course of time, 
while heavy atoms oscillate slowly, in an adiabatic hyperpotential depending on the 
proton state, through the coupling term. This separation is relevant for KHCO3 because 
adiabatic potentials for different protons states do not intersect each other. In fact, the 
separation is rigorously exact in the ground state, since protons should remain in this 
state for ever, if there is no external perturbation. Then, protons are bare fermions and 
quantum correlations arise [1]. 
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4-2. Macroscopic proton states 

Consider a crystal composed of very large numbers Na, Nf,, Nc {Af = NaNtNc) 
of unit cells labelled j,k,l, along crystal axes (a), (6), (c), respectively. Because 
centrosymmetric dimers have no permanent electric dipole, interdimer coupling terms 
are negligible [38l [39l HO] . Proton dynamics for a Ci dimer is modelled with collinear 
oscillators in three dimensions, along coordinates aijki and a2jki (a = x, y, z) with 
respect to equilibrium positions at ±aoifc«- The mass-conserving normal coordinates 
and their conjugated momentums, 

as = —1= (ai - 02 + 2ao), Psa = —7= {Pla - P2a), 

Y r (2) 

«a = ^ («! + "2), -Paa = ^ (-Pla + -^2^), 

lead to uncoupled oscillators at frequencies huga and hvaa, respectively, each with 
m = 1 amu. The difference hwsa — ^aa depends on the coupling term (say Aq,). 
The wave functions, namely '^'^jj^i{aa) / \/2^ ^n'jkii^s — "\/2«o)/"\/2 for a half-dimer, 
cannot be factored into wave functions for individual particles, so there is no local 
information available for these entangled oscillators and the degenerate ground state 
must be antisymmetrized. For this purpose, the wave functions are rewritten as linear 
combinations of those for permuted oscillators as 

-^eo,H± = 2 n ^o,H(«a) \^l,ki{(^s - V2ao) ± %,ki{(^s + V2«o)] , (3) 

where Oojfci± is the total wave function for superposed half-dimers with identical proton 
sites. Then, the antisymmetrized state vectors are: 

\Qjkl+)®\S) = |eo,H+)®^[|Tii2)-|iiT2)]; 
\Ojkl-)^\T) = \eojki-) ^ (4) 

® 77f I TiT2) + Uii2) + -^[| Tii2) + UiT2)] 

The oscillators are now entangled in position, momentum, and spin. In contrast to 
magnetic systems [H], there is no level splitting, so the symmetry- related entanglement 
is energy-free. It is also independent of the magnitude of the coupling terms between 
oscillators. Furthermore, in contrast to Keen and Lovesey [12], or Sugimoto et al. |43j . 
we argue, as an experimental fact, that there is no significant exchange integral for 
protons separated by ~ 2.2 A [H]. Protons in KHCO3 are definitely not itinerant 
particles with an energy band structure. 

The spatial periodicity of the crystal leads to extended states and nonlocal 
observables in three dimensions. Consider the vibrational wave function for a unit 
cell containing two dimer entities as Sojfcir = ©ojfcir ± ^ojkiry with r = "+" or "— ". 
Then, trial phonon waves can be written as 

, Nc Nt Na 

Sor(k) = -= 5^ 5^ XI -0^^^^ exp(zk ■ L), (5) 

1=1 fc=i j=i 
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where k is the wave vector and L = ja+ fcb + Zc, with the unit cell vectors a, b, c. For 
fermions, antisymmetrization leads to 

k ■ L = modulo 27r. (6) 

Consequently, there is no phonon (no elastic distortion) in the ground state. This 
symmetry-related "super-rigidity" [4] is independent of proton-proton interaction. 
Then, state vectors in three dimensions can be written as: 

|Ho+(k = 0))®|5); 

|So-(k = 0))®|r). ^'^ 

There is no local information available for these entangled states and the wave functions 
SoT-(k = 0) describes oscillations of the super-rigid sublattice as a whole. Finally, the 
ground state of the sublattice is: 

v^|So+(k = 0))®|^); 

VAr|Ho-(k = 0))®|T). ^ ' 

This state is intrinsically steady against decoherence. Irradiation by plane waves 
(photons or neutrons) may single out some excited states. Entanglement in position 
and momentum is preserved, while the spin-symmetry and super-rigidity are destroyed. 
However, these features are reset automatically after decay to the ground state, 
presumably on the time-scale of proton dynamics. This mechanism allows the super- 
rigid sublattice to be at thermal equilibrium with the surroundings, despite the lack of 
internal dynamics. Note that even at room temperature, the thermal population of the 
first excited proton state (7OH ^ 1000 cm~^) is negligible (< 1%). 

5. Differential cross section: calculations 

Neutrons (spin 1/2) are unique to observing the spin-symmetry of proton states ([7]). 
However, the spin-symmetry is not stabilized by any energy. It is therefore extremely 
fragile and only "noninvasive" experiments, free of measurement-induced decoherence, 
can evidence quantum entanglement [55]. For neutron scattering, this means (i) no 
energy transfer (ii) no spin-flip and (iii) particular Q- values preserving the super-rigidity. 
For planar dimers, these conditions are realized when Qx-, Qy, match a node of the 
reciprocal sublattice of protons, so that scattering events probe the super-rigid lattice 
without any perturbation of its internal state. The spin-symmetry is probed along 
the neutron-spin direction without spin-flip and the cross-section, namely the total 
nuclear cross-section an, is ~ 45 times greater than the coherent cross-section, anc, 
when the matching conditions are not realized [U HB] . Enhanced diffraction probe the 
perfect periodicity of the super-rigid sublattice, so the Debye-Waller factor is equal to 
unity. Alternatively, when the matching conditions are not realized, phonons created by 
neutrons contribute to the thermal parameters in Table [3] and to the ellipsoids in Figs 

m 

In fact, enhanced diffraction and Bragg peaks of heavy atoms are superimposed. 
Their intensities are proportional to an and o'ckco3Gxp{—2Wkco3{Q)}, respectively. 
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where o'ckcOs ~ 27.7 bn is multiplied by the Debye- Waller factor. At large Q- values, 
the contribution of heavy atoms is largely hidden by that of correlated protons, and this 
effect is more pronounced at elevated temperatures. 

The dashed lines in Fig. [1] show that proton sites are aligned along x and y, but 
not along z. Consequently, the noninvasive condition can be realized for Q^ and Qy 
exclusively, whereas Q^ never coincides with a nod of the reciprocal lattice of protons. 
Then, the diffraction pattern depends on whether we consider incoherent or coherent 
scattering along Q^. 

In order to calculate the differential scattering functions, we consider a C2/m 
super-rigid sublattice composed of primitive unit cells (a/2, b, c) with (j, k, I) indexes. 
Proton sites are at ±{xojki,yojkhO}, ^{xojkh—yojkh^}, with respect to the centre of 
split dimers. The periodicity of the grating-like structure parallel to y is -Dx/2, with 
D^ ^ a/ cos 42° ^ 20.39 A, while the periodicity of double-lines parallel to x is Dy = b. 

5.1. Super-rigid arrays in two dimensions 
The differential cross-section can be written as 



da2 ^^^ 

1=1 n T, 



dn 



Ti Tf 



E E {[exp iQy {kDy - yo) 
j=i k=i 



2 (9) 



+ TiTf exp iQy {kDy + yo)] 

X [expiQ^ (J^:e/2 - xo) + TiT/expiQ^ (j^x/2 + Xo)]Y 

X exp-2iy,(Q,), 

with A^^ = 2Na. Neutrons are scattered either in-phase {rfTi = -t-1) or anti-phase 
{TfTi = —1) by pairs of lines separated by 2xq ~ 0.6 A and 2yo ^ 2.209 A, 
respectively. The phase matching condition, namely xq (yo) commensurable with Dj;/2 
(Dy), is intrinsic to the crystal structure. {■ ■ -j^ accounts for simultaneous scattering 
by superposed singlet-like and triplet-like states ([7]), without neutron-spin flip. The 
compound Debye- Waller factor exp —21^2 (Qz), including contributions from all atoms, 
accounts for incoherent scattering along Q^. 

The diffraction pattern is composed of rods of diffuse scattering parallel to Qz, cigar- 
like shaped by the Debye- Waller factor, at Q^, Qy-values corresponding to divergences 
of (Q gathered in table [51 Such divergences occur at Qy = Uyn/yo k, Uy x. 2.86 A~^, 
since QyDy/ii ^ bUy is integer. There is no divergence for anti-phase scattering at 
Qy = ±{ny + l/2)7r/|/o, because QyDy/ii ^ 5{ny + 1/2) is not integer. For Uy even, 
QyDy/n is also even, t^ = Tf, and ridges are anticipated at Q^: = n^'^/xo ^ n^y-lO A~'^, 
since Q^D^/n ^ GSn^ is even. Alternatively, for Uy odd, QyDy/n is also odd, Xj 7^ tj, 
and ridges are anticipated at Q^ = (n^. + l/2)7r/xo, since Q^D^/n ^ GS^n^ + 1/2) is 
even. In this case, the enhanced features are best visible at A; ~ 2.57 (Fig. [3]). 
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5.2. Super-rigid lattice in three dimensions 
The differential cross-section for no spin-flip scattering events is 
da^i 



dn °^ ?? 



K Nt Nc 

Zl Zl S exp iQzlDz {[ex.piQy {kDy - yo) 
j=i k=i 1=1 

TiTf exp iQy {kDy + yo)] (^0) 



X [expiQ^ {jDx/2 - xq) + nTfexpiQ^ {jD^/2 + xq)]}^ 

Divergences occur along the previous rods of intensity at Qz = ±nz2TT/Dz, with 
Dz ^ ex cos 28° ^ 3.28 A and 2n/Dz ^ 1.92 A^^ In Fig. El these peaks clearly 
emerge from the rods of diffuse scattering depressed by the Debye- Waller factor in ([9]). 
Equations iQ and (ITO!) confirm that the enhanced features correspond to the 
reciprocal sublattice of protons represented by state vectors ([8]) based on a theoretical 
framework free of ad hoc hypothesis or parameters. The only prerequisite, namely the 
adiabatic separation, is an intrinsic property of hydrogen bonds in this crystal. By 
contrast, other models [i2l H3] cannot account for enhanced diffraction. 

6. Discussion 

Kashida and Yamamoto [2] have opposed the ordered 6*2/, to the disordered Ci structures. 
The argument consists in that for the former there is no splitting of heavy atom 
sites and the symmetrized double-wells for protons lead to quantum delocalization 
[H [HI dni [TTl HTj. By contrast, split C, dimers are supposed to be representative of 
a statistical distribution of atom positions and the double well remains asymmetric, as 
it is below Tc. 

Kashida and Yamamoto [9J or Eckold et al. [6J have observed that the P2i/a 
superlattice reflections {h + k = 2n + 1) merge above Tc into elastic diffuse scattering, 
cigar-like shaped along (a*). The intensity decreases as the temperature increases but 
the signal is still visible above 340 K. According to KY, the P2i/a local structure could 
survive above T^ as pre-transitional domains whose correlation length diverges as T goes 
down to Tc. Diffuse scattering along (a*) is invisible with SXD (Fig. [2]) but this does 
not question the existence of this signal. 

Critical phenomena have been also thoroughly investigated for the phase transition 
of KDCO3 at TcD = 353 K, taking advantage of the larger coherent scattering cross- 
section of the deuterium atom [121 [131 HI]- QENS measurements at the superlattice 
reflection (1, 2, 0) with the high resolution spin-echo technique evidence very slow 
dynamics (~ 10~^s) above Tcd and complete slowing down as T — > Tcd- These 
dynamics were tentatively associated with pre-transitional collective reorientation of 
dimers, eventually correlated with the softening of an acoustic mode propagating 
along (a*), with polarization parallel to (6) [8]. On the other hand, diffuse scattering 
around the main reflection (0, 2, 0) was attributed to the relaxation of individual dimers 

[ii[ia[ir 
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In contrast to the order- disorder schemes envisaged in previous works, the 
diffraction patterns presented in Figs [2] and [3] lead to the conclusion that the structural 
transition is of the order-order type. Then the question arises as to whether the double- 
wells undergoes symmetrization or not. 

The shape of the double-wells along the hydrogen bonds can be probed with 
Raman spectroscopy in the frequency range of the OH stretching vibrations |5l B8] . 
In Fig. [5l the submaxima observed below T^ at k, 2580, 2800, 2930, 3100 cm~^ are 
representative of the four transitions for the asymmetric double-wells of the P2i/a 
structure. Additional transitions (for example overtones and combinations) cannot be 
excluded, but a detailed assignment scheme is not necessary because the chief point of 
the argument is that symmetrization of the double-wells above T^ should be of dramatic 
consequences to the spectral profile: (i) large frequency shifts due to the collapse of 
the potential asymmetry; (ii) the symmetry-related selection rules should extinguish 
two among the four transitions. Figure [5] shows that the extended eigenstates of the 
i/OH mode probed at the Brillouin-zone-centre are identical for the P2i/a and C2/m 
structures. Symmetrization is excluded. 

The C2/m structure can be featured as follows. Firstly, enhanced diffraction 
suggests that the crystal is in a pure state {\C2/m)) composed of extended states for all 
dynamical degrees of freedom. Secondly, Raman spectra suggest that \C2/m) can be 
represented with a superposition of state vectors \Pu) and \Pd) representing P2i/a-like 
substructures, symmetric with respect to {a,c). They are labelled with "up" {u) and 
"down" (d) orientations of a pseudospin along z. According to the quantum theory of 
measurements [Tlll9l[50], the eigenstates of either |P„) or \P(i) are probed with Raman. 

For each substructure, there are two configurations, L or i?, for the proton 
sublattice. Within the framework of the adiabatic separation, the state vectors of the 
four configurations L^, Ru, L^, Rd, sketched in Fig. HI are \Pu)\HLuT)\'r), \Pu)\HRur)\T), 
\pd)\HLdT)\T), \pd)\HRdT)\T), where |p„) and \pd) represent heavy atoms. The ground 
state is a superposition of the degenerate states \Pu)\Hlut)\t) and |prf)|ifidT-)|r), while 
\Pu)\Hrut)\t) and \pd)\HiidT)\r), also degenerate, are at higher energy (AElr). The 
state vector of the crystal can be then written as 

\C2/m) = 2-y'{aiT)[\p^)\HLur) + \Pd)\HLdr)] . . 

+ (3{T)[\p^)\Hn^r) + \Pd)\HRdr)}\r), ^ ' 

with a;(T)^ + /3(T)^ = 1. The superposition of proton states is consistent with enhanced 
diffraction because the corresponding sublattices are identical in reciprocal space. 

QENS measurements of KHCO3 [6j are also consistent with (fTTj) . The inverse 
correlation time for proton transfer (tq"^ fs 2 x lO^^s"^) is temperature independent, 
without any discontinuity at Tc, because QENS, like Raman spectroscopy, probes proton 
dynamics of the P2i/a structure below T^ and the same dynamics for either |P„) or \Pd) 
above Tc. 

Furthermore, Tq^ below T^ is consistent with a superposition of proton states 
giving rise to quantum interferences whose beating frequency, z/?, ^ 2.5 x lO^^s"^, is in 
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reasonable agreement with QENS measurements [5]. Above T^, quantum interferences 
for flTTl) are cancelled by symmetry. This is at odds with the expectation that 
fluctuations should increase at elevated temperatures. However, this fluctuation-free 
state cannot be "observed" because QENS probes either |P„) or |Pd) states. 

Above Tc, diffuse scattering along (a*), very slow dynamics, and enhanced 
diffraction reported in previous works suggest dynamical inter conversion between 
superposed |P„) and \Pd) states, on the one hand, and their analogues |P2i/a)„ and 
\P2i/a)d states, on the other. This may occur through very small displacements of 
heavy atoms along (6) (Fig. [1]). Upon cooling, symmetry breaking at T^ should lead to 
either \P2i/a)u or |P2i/a)rf states, with equal probabilities. However, these structures 
are related through a translation vector (a/2,6/2,0). They are indistinguishable and 
the symmetry breaking leads to a unique P2i/a structure. 

It is often argued that a complex open system undergoing continuous interactions 
with the environment should be in a significantly mixed state that cannot be represented 
by a state vector [5T], [52], [53] . The present work leads to conclusions at variance with this 
view: the crystal of KHCO3 at thermal equilibrium with the surrounding atmosphere 
can be represented by a state vector. The gaseous environment can be treated as 
a distribution of plane waves that cannot destroy the entanglement intrinsic to the 
lattice periodicity. Only the spin-symmetry can be transitory destroyed. However, this 
reversible disentanglement does not increase with time and it is insignificant because 
the density-of-states of the surroundings is smaller than that of the crystal by many 
orders of magnitude. 

7. Conclusion 

We confirm with single-crystal neutron diffraction the C2/m structure composed of 
split Ci dimers proposed by Kashida and Yamamoto [9j, but proton positions are 
at variance with those based on X-ray diffraction. The structure of the sublattice 
of protons is practically temperature independent, apart from site occupancies. The 
enhanced diffraction pattern measured at 340 K is basically identical to that arising 
from macroscopic quantum correlations below Tc [31 [1] . 

The theoretical framework for quantum correlations is based on fundamental 
principles of quantum mechanics. The fermionic nature of protons, unveiled by the 
adiabatic separation, leads to spin-symmetry and super-rigidity. The spin-symmetry 
is probed by neutrons diffracted at Q-values corresponding to nodes of the reciprocal 
lattice of protons. The coherent scattering cross-section is increased by a factor of ~ 45, 
compared to that for uncorrelated protons, and super-rigidity cancels the shortfall of 
intensity by the Debye- Waller factor. 

Above Tc, statistical disorder, or stochastic proton jumps, or clusters, or domains, 
or uncorrelated dynamics of individual dimers, or whatever phenomena which could 
destroy the translational invariance of the crystal, are excluded. The phase transition 
occurs from one ordered P2\ja structure to another ordered C2/m structure. Both 
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phases can be treated as pure states represented by state vectors. Raman spectroscopy 
and QENS show that eigen states for protons are identical in the two structural phases. 
This suggests that the state vector \C2/m) is a superpositions of state vectors |P„) 
and \P(i) for P2i/a-like structures. This leads to the rather paradoxical conclusion 
that quantum interferences arising from superposition of proton states are cancelled by 
symmetry, as long as these states are not "observed" . 

This work emphasizes that atoms in the crystal-field of KHCO3 are not individual 
particles possessing properties on their own right. They merge into macroscopic 
states and exhibit all features of quantum mechanics: nonlocality, entanglement, spin- 
symmetry, superposition and interferences. There is every reason to suppose that similar 
quantum effects should occur in many hydrogen bonded crystals undergoing structural 
phase transitions. 
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Table 1. Neutron single crystal diffraction data and structure refinement for 
potassium hydrogen carbonate at various temperatures. A = 0.8305 A, space group 
P2i/a at 300 K (after Ref. [J) and C2/m at 340 K (this work). The dimer symmetry 
is either Ci or C2h in the two different structure models under consideration. The 
variance for the last digit is given in parentheses. 



Crystal data 


300 K 


340 K 


340 K 






Q 


C2h 


a(A) 


15.18(1) 


15.173(9) 


15.173(9) 


6(A) 


5.620(4) 


5.628(4) 


5.628(4) 


c(A) 


3.710(4) 


3.711(4) 


3.711(4) 


a 


104.67(5)° 


104.63(5)° 


104.63(5)° 


Volume 


306.2 


306.6 


306.6 


Reflections measured 


1731 


1789 


1789 


Independent reflections 


1475 


931 


931 


Reflections used 


1068 


752 


755 


ail) limit 


3.00 


3.00 


3.00 


Refinement on F 








R-factor 


0.035 


0.028 


0.038 


Weighted R-factor 


0.034 


0.019 


0.037 


Number of parameters 


56 


59 


34 


Goodness of fit 


1.025 


1.116 


1.107 


Extinction 


3260(100) 


15.4(4) 


14.9(7) 



Table 2. Atomic positions, isotropic temperature factors and site occupancies for 
KHCO3 at 340 K. The variance for the last digit is given in parentheses. 



Atom 


x/a 


y/b 


z/c 


U(iso)(A2) 


Occupation 


K(l) 


0.16526(7) 


0.00 


0.2957(3) 


0.0290 


1.00 


C(l) 


0.11951(3) 


0.5192(1) 


-0.1437(1) 


0.0186 


0.50 


0(1) 


0.19339(4) 


0.5019(1) 


0.0950(2) 


0.0354 


0.50 


0(2) 


0.0781(2) 


0.3021(1) 


-0.264(1) 


0.0337 


0.50 


0(3) 


0.0815(2) 


0.6989(2) 


-0.280(1) 


0.0301 


0.50 


H(10) 


0.0176(2) 


0.687(2) 


-0.439(1) 


0.0462 


0.25 


H(ll) 


0.0189(2) 


0.684(1) 


-0.450(1) 


0.0312 


0.25 
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Table 3. Thermal parameters in A'^ units for KHCO3 at 340 K. The variance for the 
last digit is given in parentheses. The thermal parameters account for the variation 
of the contribution of each atom to Bragg's peak intensities through the thermal 
factor r°' depending on the reciprocal lattice parameters a*, 6*, c*, and unit cell 
indexes in reciprocal space /i, fc, Z, as T"' 
2U^lha*kb* + 2U^lkb*lc* + U§llc*ha*)]. 



exp[-27r2(C/fi*/i2a 



atu2„*2 



-U^^k'^b*^ + U§iPc*^^ 



Atom 


Uii 


U22 


U33 


U23 


Ui3 


U12 


K(l) 


0.0316(4) 


0.0303(4) 


0.0227(3) 


0.000 


0.0025(3) 


-0.0003(1) 


C(l) 


0.0209(2) 


0.0117(3) 


0.0216(2) 


0.0003(2) 


0.0023(1) 


-0.0001(2) 


0(1) 


0.0247(2) 


0.0441(3) 


0.0315(3) 


-0.006(3) 


-0.0036(2) 


0.0083(1) 


0(2) 


0.029(1) 


0.023(1) 


0.041(2) 


-0.0001(8) 


-0.0067(9) 


-0.0043(7) 


0(3) 


0.030(1) 


0.028(1) 


0.0295(9) 


0.0018(7) 


0.0009(7) 


-0.0072(8) 


H(10) 


0.040(6) 


0.034(5) 


0.055(6) 


-0.003(5) 


-0.007(4) 


0.003(4) 


H(ll) 


0.034(5) 


0.023(4) 


0.036(5) 


0.004(3) 


-0.008(5) 


-0.001(3) 



Table 4. Interatomic distances in A units in KHCO3 at 340 K. The variance for the 
last digit is given in parentheses. 



C(l)-0(1) 
C(l)-0(2) 

C(l)-0(3) 

C(l)-C(l) 
0(1)~0(1) 
0(2)-0(3) 
0(2)— H(10) 



1.245(1) 
1.247(1) 
1.395(1) 
1.210(1) 
1.210(1) 
1.395(1) 
0.216(1) 
0.02(1) 
0.09(1) 
0.983(1) 



0(3)— H(10) 
0(3)-H(ll) 

H(10)— H(ll) 

H(10)— H(10) 
H(ll)-H(ll) 

0(2)-H(ll) 



1.002(1) 

1.003(1) 

0.053(3) 
0.601(1) 
0.607(1) 
0.599(1) 

0.988(1) 



Table 5. Orders n^, n^ and positions Qy, Qx of rods of intensity arising from the 
entangled array of orthogonal doubles lines of protons in two dimensions. Obs. Qx- 
positions in A~^ of the observed rods of intensity along Qz in Figs [2] and [3] QyDy/iT 
is rounded to integers. 



iy 



Qy (A- 



yDy/n k = Qy/b* Qx 



Obs. Qx (A- 






+ 1 











UxTt/xo 


0,±10 


1 


-1 


2.86 


5 


2.57 


{ux + 1/2)tt/xo 


±5, ±15 


2 


-fl 


5.71 


10 


5.14 


UxTt/xq 


0,±10 


3 


-1 


8.57 


15 


7.71 


{ux + l/2)^/a;o 


±5, ±15 
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^ • ^ 



Figure 1. Schematic view of the crystalline structure of KHCO3 at 340 K. The 
thin solid lines represent the unit cell. Top: projection onto the plane normal to c. 
The dotted lines joining protons emphasize the proton sublattice (see text). Bottom: 
projection onto {a^c). 
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Figure 3. Map of scattering intensity for KHCO3 at 340 K in between {a* ,c*) 
reciprocal planes. The arrows emphasize ridges of intensity parallel to (z) and 
perpendicular to dimer planes. 



v.v v.v 

Figure 4. Schematic representation of the Lu, Ld, and Ru, Rd, P2i/a-like 
configurations corresponding to state vectors superposed in the \C2/m) crystal state 
above Tc (see text). L and R correspond to different configurations for protons in 
dimers. Configurations u and d are symmetric with respect to {a,c). Each site is 
half-occupied. 




Macroscopic proton states 



20 



"T 1 r- 



"T 1 r- 



-T 1 r- 




2000 



2400 



1 I I r 
2800 



I I I 



3200 



3600 



Raman frequency shift (cm" ) 



Figure 5. Raman spectra of powdered KHCO3 crystal across the phase transition 
at Tc = 318 K. The spectra were recorded with a DILOR XY triple-monochromator 
equipped with an Ar+ laser emitting at 4880 A. The spectral resolution was w 2 cm~^. 
The sample was scaled in a glass tube and then loaded in a furnace. The temperature 
was controlled within ±2 K. 



